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Abstract. The eigenfunction expansion by Gegenbauer polynomials for the diffusion 
on a hypersphere is transformed into the diffusion for the Wright-Fisher model with 
a particular mutation rate. We use the Ito calculus considering stochastic differential 
equations. The expansion gives a simple interpretation of the Griffiths eigenfunction 
expansion for the Wright-Fisher model. Our representation is useful to simulate the 
Wright-Fisher model as well as Brownian motion on a hypersphere. 
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1. Introduction 

Here we show that the diffusion on a hypersphere [3j is transformed into the diffusion for 
the Wright-Fisher model with a particular mutation rate him, by using the relation, 
Xi = y 2 , where xfs denote the relative abundance of alleles and yf s denote the position 
of a particle of the diffusion on a hypersphere. Diffusion on a sphere has been applied 
to various problems in physics, chemistry, mathematics, etc. [221 [2HMU. The Wright- 
Fisher model is widely used to study the random sampling effect in population genetics 

MUM. 

We consider stochastic differential equations in the Ito sense m • We represent the 
Wright-Fisher model for k alleles as well as the diffusion on (k — l)-dimensional sphere 
ra, by using ^j) independent standard Brownian motion bij(t) i > j with the skew 
symmetry %(£) + bjft) = 0 [HI THl CHI, 18j. Our representation is useful to simulate 
the Wright-Fisher model [19, 23]. There are effective simulation methods to simulate a 
Brownian motion on a hypersphere for example [2j. Our representation is simple and 
may be useful to simulate the Brownian motion on a hypersphere. 

2. Stochastic differential equation for the Wright-Fisher model 

In the Wright-Fisher model each of the genes of the next generation is obtained by a 

random choice among the genes of the previous generation and that the whole population 

changes all at once. In Moran’s model [21] it is supposed that there are 2 N individuals 

each formed from k alleles Ai, A 2 , ..., Ak, and that each instant at which the state of the 

model may change, one individual of the alleles, chosen at random, dies and is replaced 

k . 

by a new individual which is Ai with probability 7^7 j w here ki is the abundance of 
the allele At. It is supposed that the probability of any individual ’’dying” during an 
interval (£, t + dt) and then being replaced by a new individual is Xdt. Hence the mean 
number of such events in unit time is A and the mean length of a generation is A -1 . The 
following model of interacting particles [T2] is equivalent with the Moran model. 

Consider a population of N particles each of which is one of k types, Ai, A 2 ,..., A^. 
The types may represent species, alleles, genotypes or other classification. We then 
consider random interactions between particles, which are assumed to occur at the rate 
Xdt per time interval (t, t + dt) for each particle. If a pair of particles of different types i 
and j interact, then after the interaction the both particles are the type i with probability 
1/2 and the type j with probability 1/2. If the type of the interacting particles are the 
same, no change occurs. 

We can approximate the behavior of our interacting particle system by a stochastic 
differential equation (JT]) . I 11 it, the relative abundance of type i increases by 

Cy/xi(t)xj(t ) dbij(t) and decreases by cyjxi(t)xj(t) dbijft ) with the interaction of the 
particles of type j, where c = ^JX/2N. Hence our interacting particle system 
automatically makes the following equation ([I]), which has the genetic drift matrix with 
the elements c 2 Xi(t)(Sij — Xj(t))dt for i,j = 1, 2,..., k as covariances. 





Diffusion on a Hypersphere: Application to the Wright-Fisher model 


3 


For i,j = 1,2consider 



( 1 ) 


3 = 1,*¥=3 


with bij(t) + bji(t) = 0, where bij(t) (i > j ) are mutually independent one-dimensional 
Brownian motion with the mean 0 and the variance t nn Eina ini- This representation 
of the Wright-Fisher model mi is applied to make computer simulations on models in 
population genetics [IHl 1233 • Here we make use of this equation ([!]) to discuss the exact 
solutions of the diffusion equations considering the diffusion on a unit sphere. 

3. Diffusion on hypersphere 

The isotropic diffusion on the (k — l)-dimensional unit sphere is given by the following 
stochastic differential equation in the fto sense. 


dyiit) = -£-(k - 1 )yi(t)dt + ^^2yj(t)dbij(t) 


( 2 ) 


with bij(t) + bji(t) = 0, for i = 1, 2,..., k, where bij{t) (i > j ) are mutually independent 
one-dimensional Brownian motion with the mean 0 and the variance t. Let dy = 
(dyi ,..., d/ijk). By using the fto calculus [20], we have for the stochastic differential 


d(5>*(*) 2 ) = o. 


(3) 


Hence starting from a point on the unit hypersphere, the trajectory is on the unit 
hypersphere J2iUi(t) 2 = 1. Consider the dot product for a unit vector 1 = (Zi, Z 2 , ...,4), 


2 

1 • dy = E l & E(^(*) d M*)) + ~4~( k - 1 )Vi{t)dt\. 

i j 


(4) 


On the tangent hyperplane at y = (jq, y- 2 ,..., yk) of the hypersphere, for 1 which is 
perpendicular to y, we have 


ndy = | E ME %■(*)<%(*))■ 


i j 


(5) 


We have the expectation 

£(l-dy) = 0. 


( 6 ) 


We have 


Var{ 1 • dy) = (^) 2 Ete ~ hyi) 2dt = 


(7) 


i>j 


since 


Eto ~ Vih ) 2 = Et(M/j ) 2 + (ydj )' 2 ~ zhyjyilj] + (EM /*) 2 



= Yiihyj) 2 + C Vih ) 2 ] + Y&yi) 2 = (E*?)(E tf) = 1- 



( 8 ) 
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Var (1-dy) does not depend on the direction 1 at the tangential point y of the tangential 
plane. Hence we see the diffusion is isotropic and we have the Fokker-Planck equation 
on the hypersphere for the Laplace-Beltrami operator [3] in space S k - 1 , 

y, t\y', 0) = DA Sk ~ 1 p(y, t\y, 0) (9) 


where D = 
<9i, 6> 2 , - - -, Ok -1 


8 


from (JTJ). The operator A Sk ~ 1 is for k — 1 spherical coordinates 
of the unit vector y = (yi,...,y k ) of S k -i dehned by the relations 


yi — sin Ok-i-.. sin0 2 sinffi (10) 

y 2 — sin Ok-i-.. sin 0 2 cos 0 i (11) 

: ( 12 ) 

Vk -1 = sin Ok -1 cos Ok- 2 (13) 

y k — cos 0 k -\, (14) 


where 0 < 0i < 2 tt and 0 < 0i < n for i ^ 1. The integration measure in S k - 1 is dehned 
as dy = -r— sin^ -1 0 k -\... sin 0id0i...d0 k -i , where A k -\ = 2 n k ^ 2 /T(k/ 2 ) is the surface 
of the sphere S k -i- With this normalization J s dy = 1. 

The solution [3], the transition probability density for the isotropic diffusion on 
hypersphere, is obtained by using the definitions and notations of the chapter IX of 
the book [24]. The Gegenbauer polynomial C P L is dehned as the coefficient of h L in the 
power-series expansion of the function 

OO 

(1-2 th + h 2 )~ p = J2 C P L (t)h L ■ (15) 

L =0 

The transition probability density from y' (with n — l spherical coordinates Of ..., 6 ) ( i _ 1 ) 
at 0 to y (with 0 1 ,..., Q n -\ ) at t > 0 is given [3] by 


p(y,£|y',o) 


£ 2 L ~k- 2 2 ^ /2 ~ ex P (~ DL ( L + k- 2 )t), 


A-k-i L=0 

for the dot product y • y ' — UiUi , where 


dWh) = E (-in 

3=0 


LL/2J r(L-j + k/ 2 - 1 ) 


T(k/2-l)j\(L-2j)\ 


ET( 2 -) 


L-2j 


(16) 


(17) 


4. Solution for the Wright-Fisher model with parent-independent mutation 

Let A 1 ,...,A k denote k allele types in a population. The general neutral alleles model 
has a probability ii y - of a mutation from Ai to Aj. An exact solution is known for the 
case of the 0 mutation rate case in three dimensions [16]. Assuming parent-independent 
mutation [TJ , Uij = Uj. an expansion of the transition probability density [7], from the 
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initial relative abudances x' = [x\,..., x' k ) to the relative abundances x = (aq, ...,Xk), is 
obtained for the following Fokker-Planck equation (fT8]h 


<9p(x, t | x', 0) 
dt 

where Mj = — px i: p = 

|23 [25] is given by 


h dx i 2 h 9x j 

E|Li q- For > 0,i = 


— - ^Mi)p(x,t | x',0) (18) 

1,2the stationary density 


rMn*r'm («) 

1=1 

The solution of equation (fT8|) is given by an expansion in orthogonal polynomials [7]. 
For the case e* = e, i — 1, k, we have the solution [7], 


P(x,t | x',0) 


k £— 1 oo 

kL A 


= r (» II Six (E exp{-Jn(n - l)t - Erai}Q n (x, x')) 


»=1 ^E) n=0 

for Xj > 0, j = 1, 2,..., k, where Q n is given by 


( 20 ) 


Q n (x,x') = (p + 2n- l)(n!) 1 (-1)’ 

m =0 


n 


. m 


(P T Tn)(n— l)6n (21) 


with 


6n /^(m)F(£) 


E 


m. \ x j x j 


,/ \L 


-hr - 

■ /■—■ i 


( 22 ) 


/i+...+Jfc=m,o<i 3 -,j=i,...,fc W 3 F(Zj+e) 
for the notation a?( m ) = a(a + l)---(a + m — 1). We give the particular cases as, 

Qo(x,x') = l, (23) 

Qi(x,x') = (p +1)(6 - 1), (24) 

(25) 


f? 2 (x, x') — ~(p + 3)[(/i + 2)6 — 2 (p + 1)6 + //)], 


with 


6 = h-E 


x j x j 


3 = 1 

k f ry* _ sytf \ 2 2yi _ /y» , 2V>^ 

6-^ + i)E(-T^4 + 2 E 


i=l e ( £ + 1 ) 




(26) 


(27) 


For i = 1,..., k, putting y\ = x i: the isotropic diffusion ((2D on (k — l)-dimensional sphere 
is transformed to 

dxi{t) = — (1 — kxi(t))dt + ^2 cJxi(t)xj(t)dbij(t). 


j= 1 


The diffusion process governed by equation (|T8l) is expressed as, 


( 28 ) 
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Thus the corresponding Fokker-Planck equation of Eq. (l251) is a special case of the 
equation (fTST) . putting £, = 1/2 for all i and taking c = 1. 

The transition probability density p(x, t | x',0) for equation ([T8l) with £i = 1/2, 
i k on the relative abundances x = (aq,..., x fc ) starting from x' = (x' 1 ....,x' k ) is 

represented by using the transition probability density p(y,t | y',0), for equation (fT6l) . 
on y = with spherical coordinates 9,,...,9 k of (k — l)-dimensional sphere 

starting from y' = ( y[ , ...,y' k ) with spherical coordinates 9(, ...,9' k , as 


p(x, t | x', 0 ) 


<9(oq,. 

“5 %k — l) 

d(9 u . 

• Ok- 1 ) 


E 


'2 _ „/ n ^„ ( / „, 2 _ 


We have 


d(xi ,. 

• 5 3'n— l) 

d(9 u . 

■ 9 n - 1 ) 


= 2 


d9,...d9 k _, 


p{y,t I y',o) 


d(yii -1 Vk-l) 



Uk- 1 ) 

0(9,,. 



d9\...d9 k -\ 


i=l 




(30) 

(31) 

(32) 


Considering equation iITtTli (??), dTTll (IH71) we have 


p(x,( | x', 0 ) 

r « :/2 >i]>r 1/2 ) 2 - 


7T 


k/2 


E 


i=l 


:2 _„/ n o,2_ 


00 9 T 4- h — 9 / 

E —r—o— c L 2 ~ 1( (y • y') exp (~dl(l + k- 2 )t). 


L =0 


k — 2 


(33) 


Thus we proved the solution (133|) is equivalent with the solution (!20|) for the case 
s = 1/2. Let us observe the two forms of the solution. Consider the eigen values of the 
two diffusion equations (|9]) and (|T8j) . The Caillol solution [3] of (EJ), expanded by the 
Gegenbauer polynomials, are applied in equation (fTTI) . The Gegenbauer polynomials for 
odd L in equation (133]1 cancel with each other as we can see from equation (fI71) . Hence 
there is no term for odd L in equation mm- Put L = 2n, and take D = 1/8 (c = 1 in 
(ED), for particular mutation rate £ t = 1/2, i — 1,..., k. Then the exponential functions 
exp(— DL(L + k — 2)t) in equation (133]) coincide with those of Griffiths solution ([20]) . 
The above <5 o( x j x/ )) Qi( x > x ')> an d Q 2 ( x , x/ ) hr the solution ([20]) are obtained from the 
terms in the expansion ([33]) by using ay = yf and x\ = y'f for i = 1,..., k, as we can 
calculate easily. Also we see 

k „ k rr £ ~ 1 

(34) 


r(i/2) ;nA: 


7 r 


k/2 


k T .e-1 

= r(M)n 


i .1 r(e) 

Our study naturally connects the diffusion on a hypersphere with the Wright-Fisher 
model for a particular mutation rate of the parent-independent mutation. 

Remark 1. The connection between diffusion on a hypersphere and the 1- 
dimensional Wright-Fisher diffusion has been derived before (p338 in [15] ). The 
equation (13.33) in [15] corresponds to equation (TT 6 ]) . They derive the spectral expansion 


















Diffusion on a Hypersphere: Application to the Wright-Fisher model 


7 


for a wider range of parameters. To extend this argument to the multidimensional case 
is one of our next problems. 

Remark 2. The eigenfunctions in our equation (T33T) and the eigen function (T2Tj) in 
the Griffiths solution for the particular mutation rate look different, although they are 
identical. We do not have the algebraic proof in our present work and wish to have it 
in our next work. The form of the polynomials in our expansion means that they are 
the reproducing kernel polynomials for orthogonal polynomials on the Dirichlet measure 
[TO] , because all pairs of polynomials of the same degree are added to form the coefficient 
of a given eigenvalue . 
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